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ON THE IRREVERSIBLE DYNAMICS EMERGING 
FROM QUANTUM RESONANCES 

M. KONENBERG AND M. MERKLI 


Abstract. We consider the dynamics of quantum systems which possess stationary states 
as well as slowly decaying, metastable states arising from the perturbation of bound states. 
We give a decomposition of the propagator into a sum of a stationary part, one exponentially 
decaying in time and a polynomially decaying remainder. The exponential decay rates and the 
directions of decay in Hilbert space are determined, respectively, by complex resonance energies 
and by projections onto resonance states. Our approach is based on an elementary application 
of the Feshbach map. It is applicable to open quantum systems and to situations where spectral 
deformation theory fails. We derive a detailed description of the dynamics of the spin-boson 
model at arbitrary coupling strength. 


1. Introduction and main result 

1.1. General setup. Let Lq be a self-adjoint operator on a Hilbert space "H and consider 
(1.1) L = Lq AJ, 

where A G R is a perturbation parameter and / is a self-adjoint operator. It is assumed that L 
is self-adjoint. We suppose that the spectrum of Lq is absolutely continuous (possibly, but not 
necessarily semi-bounded) and that Lq has hnitely many eigenvalues e with hnite multiplicities 
nie- All the eigenvalues of Lq are embedded in the continuous spectrum. The general problem 
we consider is how the stability or partial stability or instability of these embedded eigenvalues 
under the perturbation affect the dynamics generated by L. One may readily incorporate into 
our results and proofs the case where Lq has also isolated eigenvalues by using ordinary analytic 
perturbation theory on the corresponding subspaces. 

In the setting of usual analytic perturbation theory |23] an isolated eigenvalue e of Lq with 
multiplicity rUe splits, under perturbation, into a group of eigenvalues ..., of L (1 ^ 
ie ^ LUg), in the sense that Ef.j = Eej(A) —)■ e as A —)■ 0, for j = 1,... ,if.. For hxed e, the 
sum of the multiplicities of the eigenvalues Ef,j equals mg. On the other hand, it is well known 
that embedded eigenvalues can be unstable, or partially stable, under perturbation. Instability 
means that L does not have any eigenvalues in a neighbourhood of e for small A. Partial 
stability means that the embedded eigenvalue e of Lq splits into a group of eigenvalues of L 
whose sum of multiplicities are strictly smaller than that of e. 

In this paper, we consider the situation where Lq has unstable and partially stable eigenvalues, 
and where the partially stable ones undergo a reduction to dimension one under perturbation. 
Namely, close to any eigenvalue e of To, the operator L either does not have any eigenvalue (e 
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unstable) or L has exactly one simple eigenvalue E^, close to e, meaning that limA^.o-^'e = e. 
It is supposed that all eigenvalues of L are of this form. One may develop the arguments of 
the present paper in the more general setting where close to every e, L has several eigenvalues 
El,... ,E(^ and each of them may be degenerate. We do not do this here to keep the exposition 
simpler. 

The dependence of Ee on A is not governed by usual analytic perturbation theory, since the 
unperturbed e is an embedded eigenvalue of Lq. However, suitably modihed expressions from 
analytic perturbation theory of isolated eigenvalues will still play a role in the present setting. 
Let Pg be the spectral projection of Lq associated to the eigenvalue e. If e was an isolated 
eigenvalue of Lq, then the hrst and the second order corrections (in A) of eigenvalues would 
be given, according to analytic perturbation theory [23], by the eigenvalues of Pg/Pg and of 
Pg/Pg-^(Lo — e)“^/Pg, respectively. 

We assume that 

(Al) For all eigenvalues e of Lq, 

(1.2) PgJPg = 0. 

For embedded e, the resolvent P^{LQ—e)~^ does not exist as a bounded operator, so PgJPg*-(Lo~ 
e)~^IPe is not dehned, typically. Nevertheless, we can replace e by e — ie and consider e small. 
This suggests that the second order eigenvalue corrections to e are linked to the level shift 
operator 

(1.3) Ag = -PgJPg^(Po - e + iO+)-'JPg, 

where i0+ indicates the limit of the resolvent {Lq — e + ie)“^, as e —?■ 0+. The existence of 
the limit is guaranteed by assumption (A2) below (take A = 0 in the resolvent in fll.Sp i. The 
operator Ag is represented by an mg x mg matrix. 

Let Q be an orthogonal projection and denote 

(1.4) R. = {L- z)-^ and Pf = {Q^LQ^ - z)-^ (RanQX . 

In the following, we denote by C{(f, fj) a constant which is independent of x and A but which 
may depend on 0, '0 G "H. 

(A2) (Limiting Absorption Principle.) There is a dense set P <Z Li with RanJPg C P (Ve) 
and there is an a > 0 such that the following hold. 

(1) Let Se = {z E C- : |Rez — e| ^ a}. Here, C_ denotes the (open) lower complex 
half plane. For all e and all eP, we have 

(1.5) sup ^ C'(0,'0) < oo, k = 0,...3. 

zeSe 

(2) Let Poo = {2; e C_ : |Rez — e| > a for all e}. For all eP, we have 

(1.6) sup |^(0,P^V’)l ^ C'(0,'0) < 00, k = 0 ,l. 

ZSSoo 

Assumption (A2)(l) implies that Pf.IR^^'IPe extends as a twice continuously differentiable 
function to real z with \z — e\ ^ a. Moreover, the estimates fll.5p and fll.6l) with fc = 0 imply 
that Pf-LPf- and L have purely absolutely continuous spectrum in the interval {e — a,e + a) 
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and in the region {x G R : |a: — e| > a for all e}, respectively (see e.g. PI or [lO], Proposition 

4.1). 

(A3) The operators -^PelIPg is bounded, uniformly for 2 ; G C_ with |Re 2 ; — e| ^ a and 
for all |A| ^ Aq. 

The Feshbach map associated to an orthogonal projection Q, applied to L — z, is dehned by 

(1.7) d{L-z-,Q) = Q{L-z-LRfL)Q. 

It follows from assumptions (Al) and (A2) that 

hm d{L -Eg + ie; Pg) = d{L - Eg- Pg) = Pg{e -Eg- A^IR^l I)Pg. 

Let V'Ee = Le s.t. = Pg-ipE^- Note that Pg-ipE^ 7 ^ 0 for otherwise P^LP^-ipE^ = 

EgPj-'ijjEe, which cannot hold for small A, since P^LPj- has purely absolutely continuous 
spectrum in a neighbourhood of e due to Condition (A2)(l). We normalize -ipEe as UPeV’Cell = 1- 
By the isospectrality property of the Feshbach map (see Appendix 0 Proposition IB.2p . we 
have 

(1.8) 0 = diL- Eg-, Pg)Pg^E. = (^ - PgIRlUo^IPg)Pg^E.. 

This, together with Conditions (A2)(l) and (A3), implies that .^a := is bounded in A 
for small A and that PglPg = Ag + 0(|A| + |e — Eg\) = Kg + 0(|A|). On a suitable 
sequence A„ —)■ 0, we have ^a^ ^0 and Pg'ipEX^n) PeV'o for some ^0 ^ R and some unit 
vector Pgijjo (Bolzano-Weierstrass). Consequently, taking A ^ 0 in (II.8p along this sequence 
gives 

(1.9) AePeV'O = ^oL’eV'O, 

showing that ^0 is a real eigenvalue of Ag. For ease of presentation, we assume the following. 

(A4) (Fermi Golden Rule Condition.) The eigenvalues of all the level shift operators Ag are 
simple. Moreover, 

(1) If e is an unstable eigenvalue, then all the eigenvalues Ae,o, • • •, K,me-i of Ag have 
strictly positive imaginary part. 

(2) If e is a partially stable eigenvalue, then Ag has a single real eigenvalue Ae,o. All 
other eigenvalues Ae,i,..., Xg,me-i have strictly positive imaginary part. 

Under condition (A4)(2), the set {^a = for A small has a unique limit point ^0 and 

we have Ae,o = ^o- Having only simple eigenvalues, Ag is diagonalizable and has the spectral 
representation 

me —1 

(1.10) Ae = Ae,,Pe,„ 

j=0 

where Pgj are the (rank one) spectral projections. We introduce the notation 

(1.11) a -<b, 
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where a is a complex number, a vector or a bounded operator and 6 > 0, to mean that 
|a| ^ const.6, where | ■ | is the appropriate norm and const, is a constant which does not depend 
on the coupling parameter A, nor on time t. 

Theorem 1.1 (Resonance expansion of propagator). There is a constant c > 0 s.t. for 0 < 
|A| < c the following holds. Denote the projection onto the eigenvalue of L by IIe^. Let 
t > 0, (fj'ip ^ D s.t. L(j), L'0 G D. Then 

me —1 

( 1 . 12 ) = Y. 

e partially stable j=l 

me —I 

+ Y Z <*‘“'^^’“•’'’('>■"1,^)+«(*)■ 

e unstable j=0 

where 

(1.13) R(t) -< 1. 

The exponents Oej and the operators - are close to the spectral data of the level shift operator 
Ae, (ll.lOp . Namely, 

(1.14) a.,, = Kj + 0(A), n),, = P.J + 0(A). 

Remarks and discussion. 

1. For an expansion of for f > 0, simply take the adjoint of (I1.12p . 

2. The exponents a^j are the eigenvalues of an explicit matrix. They can be calculated to 
all orders in A (see Lemma l3.ip . The operators flg^- have also expressions calculable to all 
orders in A (see flS.Sp l. Those “complex energies” are called resonances. They coincide 
with the eigenvalues of the spectrally deformed generator of dynamics in situations where the 
latter exists. 

3. The remainder term is small relative to the contributions of the exponentially decaying 

terms in fll.l2p for times t satisfying ^ 1/f, where 

(1.15) 7 = min | ImOe j }. 

The inequality C ft, for some (large) C is equivalent to ^ For small A, 

it is valid for intermediate times, to < t < ti, with to = C + 0(7A^) and ti ~ 1/(7A^). During 
this time-interval, the decay of fll.l2p behaves as exponential, to leading order. 

4. The Fermi Golden Rule Condition (A4) guarantees that instability (and partial stability) 
of eigenvalues is visible at lowest order, O(A^), in the perturbation. It may happen that 
resonances acquire non-vanishing imaginary parts only at higher orders in A. Our method can 
be adapted to describe this situation. The dynamical consequence is a slower decay of the 
corresponding directions in Hilbert space, see also [3T] . 

5. Increasing the regularity assumptions on the vectors 0, 0 allows to show a faster poly¬ 
nomial decay of the remainder than fll.131) (this amounts to taking higher z derivatives of the 
resolvent, c.f. fll.lTp l. 















ON THE IRREVERSIBLE DYNAMICS EMERGING 


FROM QUANTUM RESONANGES 


5 


1.2. History, relation to other work. The analysis of resonance phenomena has a long his¬ 
tory and plays an important role in quantum physics p El ESI an . Its modern description, 
involving dilation analytic Hamiltonians ([HE]), was given in [38] and further developed in 
[39l [371 [20] . We refer to [19] for a textbook presentation and many more references. A time- 
dependent theory of quantum resonances was established in [32], inspired by [iQ] and further 
developed in [25]. In these works, as well as in m, a variant of the Mourre theory in combi¬ 
nation with the Feshbach projection method is used to link dynamical properties of quantum 
systems to spectral objects. The approach of the present work is, in spirit, similar to m-, 
see the end of this paragraph for a comparison. Nevertheless, all the above-mentioned works 
require regularity conditions within Mourre theory that do not allow the treatment of open 
quantum systems at positive temperature. In the context of open quantum systems at posi¬ 
tive temperature, the link between quantum resonances and approach of an equilibrium state 
has been pioneered, using complex deformation theory, in [^ [7] . The work [7] is based on a 
sophisticated renormalization group method initiated in [HUE]. Recently, a method based 
on graph expansions of the propagator rather than purely spectral considerations was given in 
[T2] . The spectral approach has been further developed to yield a detailed description of open 
systems dynamics in terms of resonances in [331 EH ES] , with applications to quantum informa¬ 
tion theory [231 ED] and quantum chemistry [22] . The spectral analysis and its consequences for 
“return to equilibrium” based on Mourre theory and positive commutators was carried out in 
[131 EH EH Hi]. However, these papers are limited to the study of the spectrum of the Liouville 
operator with the goal (typically) of showing that it has a single, simple eigenvalue at zero (and 
absolutely continuous spectrum otherwise). This information alone does not provide any detail 
about the dynamics other than ergodicity. However, one is interested in information such as 
directions of decay and decay rates which describe, for example the speed of thermalization, 
decoherence and the dynamics of entanglement. In the method of complex deformation, com¬ 
plex resonance energies are linked “automatically” to the decay rates of reduced density matrix 
elements [31] • The same expressions describing those decay rates appear as well in Mourre 
theory as a consequence of the Fermi Golden Rule (see also d). however, linking them to the 
dynamics, and in particular to time decay, is more delicate and has not been done previously 
for open systems. We show in the present paper how to extract the detailed dynamical infor¬ 
mation from the Mourre theory in a technical setup that includes positive temperature open 
quantum systems. It is important to be able to handle these questions using a softer approach 
than the spectral deformation one. Indeed, the applicability of the latter demands much more 
regularity from the models and, in some physically relevant situations, the spectral deformation 
technique is not applicable at all. This happens for the spin-boson model at arbitrary coupling, 
whose ergodicity has been shown recently in [21] using Mourre theory. As an application of our 
method, we give a detailed expansion of the propagator of this model in the present paper. 

The philosophies of [11] and the present paper are similar, in that the common main idea is to 
write the propagator as a contour integral over the resolvent and subsequently use the Feshbach 
map to analyze the latter. However, right from the start, the technical assumptions are very 
different. A core assumption of m is that multicommutators of the Hamiltonian H with the 
(Mourre theory) conjugate operator A are relatively FT-bounded. While this is typically true 
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for, say, for Schrodinger operators, it is not so for open quantum systems. The problem comes 
from the fact that the number operator is not bounded relatively to the free held Hamiltonian 
(so already relation (2) of [11] is not valid). The situation even is worse for the spin-boson 
model at arbitrary coupling, where each successive commutation of the Liouvillian with the 
conjugate operator produces a more singular operator, as explained in [ 2 l| (this is the reason 
why the spectral deformation theory fails). 

The dissipative character of the system is guaranteed in [11] by assuming a Mourre estimate, 
localized spectrally on a subspace around the embedded eigenvalue in question. Accordingly, 
the main result of im describes the dynamics of an initial state (wave function) which is 
spectrally localized close to the embedded eigenvalue. This makes good physical sense in 
the context of, say, Schrodinger operators, where long lived initial states are expected to he 
close to unperturbed bound states. However, in open systems problems, one considers initial 
states which are spatial perturbations of equilibrium states and which are not at all spectrally 
localized relative to the Liouville operator. Our dissipation assumption (A2) is thus a Limiting 
Absorption Principle which is not spectrally localized and which produces results for initial 
vectors which are not spectrally localized. Condition (A 2 ) can be heuristically understood in 
the context of open systems as saying that the reservoir stays essentially in its equilibirum state 
during the dynamical process (the Born approximation). 

We mention that the main result of HH is stated and proven for a simple unperturbed eigen¬ 
value. In open systems however, the origin is always a degenerate eigenvalue of the unperturbed 
Liouville operator and so we have put in place a formalism that works for the degenerate case 
as well. 


1.3. Outline of the proof of Theorem II.11 For ^|J G dom(L), we have ([IS], Corollary H 
3.6) 


(1.16) 


= - 

^ 27ri 


dz, 


R—iw 


where tc > 0 is arbitrary and Rz = {L — z) \ see (II. dh . We write e'*^ = and integrate 

by parts in fll.l 6 p to obtain 


(1.17) 


1 1 


^ it 27ri 


The boundary terms in the integration by parts vanish since Rzip 0 as |z| —?■ oo. 

We analyze separately the contributions to the integral in fll.lTp coming from in different 
regions on the line of integration. Dehne the gap of all the clusters of resonances by 

(1.18) (5 = min {|Ae,i - Aejj : i ^ j} > 0 




and denote the eigenvalue gap of Lq by 

(1.19) = min{|e — e'l : e 7 ^ e'} > 0. 

e,e' 

Let 


( 1 . 20 ) 


0 = 1 mm 


{c5, g, ImAej : (e, j) s.t. ImAej > O} > 0. 
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Here, c is a constant not depending on A,w (its origin is explained in Lemma Id.ljl . For any 
eigenvalue e of Lq, set 


(1.21) Qe = {x — iw : |x — e| ^ a} 
and set 

(1.22) goo = {x-iw : a:eR}\Ue^e. 
It follows from fll.171) that 


(1.23) {(j), ^ Jeit) + Joo(t), 

e 

where 

(1.24) = iR.i,}dz. 

We now apply a suitable Feshbach map to the resolvent Rz in (ll.24p . with a projection depend¬ 
ing on the region of integration. Let P be an orthogonal projection and recall the notation 
fll.dp . The resolvent has the representation 

(1.25) Rz = diz)-^ + B{z) + R^, 
where ( 5 ^( 2 ;) = ^{L — z,Q), see fll.7p . and 

(1.26) B{z) = -d{z)-^QLRf - RfLQdiz)-^ + RfLQd{z)-^QLRf. 

We explain these relations and some properties of the Feshbach map in Appendix [Bl For z ^ Qe, 
we choose the projection Q in the Feshbach map to be Pg- For G Qoo, the argument is simpler, 
see Section [331 

Let us assume that the unperturbed, partially stable eigenvalue of Lq is at the origin, e = 0. 
(Otherwise see section (21) Then L has a simple eigenvalue E = Eq with P —)■ 0 as A —)■ 0. To 
analyze Jo(t), we write, according to fll.25p . 

(1.27) Jo{t) = ^ {((^, f 5^(2)"V) + (0, £B{z)iIj) + (0, £R^°'ip)}dz. 

The Feshbach term is ^^(. 2 ) = —z + where Az = —PqIR^°IPq. For z = 0 and A = 0, 

Az is just the level shift operator Aq, (11.311 . We show in Lemma [3.11 that Az is diagonalizable, 
Az = ^o,j{z)Qj{z), and that the eigenvalues aQj{z) of Az satisfy ao,o(P) = P/A^ for all 

A 7 ^ 0 (this follows from the isospectrality property of the Feshbach map and the fact that P 
is an eigenvalue of L) and aoj(^) = ^o,j + 0(A^ + \z\), j = 1,... ,me — I (since Az is close to 
Aq). Then we can write 


A 

dz 




-1 



j=0 


Qjjz) 

-z + A‘^aoj{z)' 


(1.28) 
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We are interested in the singularities of this function as is close to the real axis. They come 
from the denominator. To understand the nature of the singularities, and since 2; 1—)■ Qj{z) is 
regular, consider Qj{z) ~ Qj(0) for a moment. Then 

me— 1 


(1.29) 


j=0 




^ (z-A2aoj(^))' 


■Q,(0). 


For j = 0 we have —E + A^ao,o(-E') = 0 (see above) and the corresponding summand is 


1 


1 - A^ 


0,0 


.z 


{z - Ey {1 - A2[ao,o(2:) - aofl{E)]/{z - E)y 

1 1 - 


{z - Ef (1 - A^a',^,{E)y 


Qo(o) 


Qo(o) = 


Qo(o) 


{z - Ey 1 - A^E,^,{Ey 


By using that the projection associated to the eigenvalue £' of L is given by 

(1.30) n^; = lim (ie)(L — i? + ie)“^ 

£^>0+ 

and decomposing the resolvent in this limit according to fll.25p with projection Pq, we identify 

(see (1328) 

Q«(E) „ „ p 

1 -AV„,„(E) ” ”■ 

For j > 0 we have aoj(O) = Xqj + O(A^) which is in the open upper complex half plane and 


the corresponding summand in (11.291) is 

1 - 


;QM 


Qj{0) 


{z - A‘^ao,j{z)y (1 - A2aoj(0))2 

In Section [3] we make these arguments rigorous. Namely, we show that 

1 _Q,(0) 


(1.31) 


=(,_b,s 


PallEPo + T— 


^ (^-A2aoj(0))^ 


+nz), 


where e^^^T{z)dz -< 1. Now we have to multiply fll.dip by e‘*^ and integrate over z E Oo = 
[—a, a] — iw. Having in mind a standard argument from complex analysis, we complete the 
path Qq into a closed contour (a rectangle with a ‘roof’ parallel to Qo but shifted far into the 
upper complex half plane). We then use the Cauchy formula for contour integrals to get 

/* 


(1.32) 


= 27ri ■ ife‘* +0{l/t) 


Jgo - Ef 

The 0{l/t) term is the contribution from the parallel vertical sides of the rectangular closed 
integration path (see (I3.2ip ). In a similar way, we treat the sum in (11.311) . Here the poles are 
at z = A^aoj(O) and so 

f 

(1.33) 


'go (^-^^«o,j(0)) 


;dz = 27ii-ity^^ +0(l/f). 
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Combining (ll.31j) with (ll.32p and (ll.33p yields 

^ me-l 

(1.34) ^ e'‘’^{4,,P„YlEPai:)+V 

Next we deal with the second integrand in fll.27p . Using again the spectral representation of 
5(2 ;)we have from fll.26p 


(1.35) 





j=0 


-z + A‘^aoj(z) ’ 


where 


(1.36) q,(z) = -A [Q^(z)PoIRp + R^°IPoQ,iz) - AR^°IPoQ,{z)PoIR^<^] . 

The expression fll.351) has the same structure as fll.281) . We readily obtain, in analogy with 

(Ol, 

H H me-l 


i=i 


Proceeding as above, after (ll.30p . we identify PoIfETb + [e) — (see also (I3.43p ). 

Finally, since by Assumption (A2)(l), 


JGo 

we obtain 

me —1 

(1.37) Mt) = {,p,iIe4>) + 5^e'“’‘».ji“>(.).,ni,y) + o(iA), 

i=i 

where IIq^ = (5o(0) + 0(1^1) = Po,j + 0(1^1)- This explains the contribution of a term on the 
right side of fll.l2p coming from a partially stable eigenvalue e (= 0). The analysis for unstable 
e follows using the same arguments. Finally, to deal with Joo{t), we write 

(1.38) RziIj = {z + i)~‘^Rz{L + i)^i/i - {z + i)“V - (^ + i)~^(-^ + 

which is valid for ijj G dom(L^). The negative powers of .2 help the convergence of the z-integral 
over Qao- The bound Joo(^) ~< lA is then easily reached using (A2)(2). 


2. Application to open quantum systems 
2.1. Setup. The Hilbert space is the product of a system and a reservoir part, 

( 2 . 1 ) H = Hs^Hk. 

The self-adjoint generator of dynamics, called Liouvillean, is of the form fll.ll) . where Lq, the 
free (non interacting) Liouvillean, is a sum of a system and a reservoir contribution, 

( 2 . 2 ) 


To — Ts -I- Tr, 
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and / is the system-reservoir interaction operator. We consider the system to be hnite- 
dimensional and the reservoir to be an inhnitely extended free Bose gas at positive temperature, 
as we explain now. 

Let S be a quantum system with pure state space i^s of dimension do < oo. For instance, for 
a spin 1/2, do = 2. Then the Hilbert space "Hs in 112.ip is the GNS space (Liouville space) 

(2.3) •Hs = -^s®%, 

so that d = dim'Hs = d^. The doubling of the pures state system Hilbert space in fl2.3p 
allows to represent any (pure or mixed) state of S by a vector. Namely, let p be a density 
matrix on If has the diagonalized form p = f^ which we associate the vector 

'hp = Ei ® i/j e ils G) i^s (complex conjugation in any hxed basis - we will choose the 

eigenbasis of the system Hamiltonian). Then Tr(pH) = (Tp, (H ® ls)fllp) for all A G B{S^s) and 
where Is is the identity in i^s- This is the GNS representation of the state given by p [H IM] . 
Let Hs = Ej d^j\Pj){Pj \ ho the Hamiltonian of S, acting on i^s- The equilibrium density matrix 
is ps = e“^'^®/Tr which is represented on Hs by the vector 

(2.4) Gs,/ 3 = (Tre“^-^®)“^ (g) ipj. 

3 

The (GNS) Hilbert space of the spatially inhnitely extended free bose gas, for states normal 
w.r.t. the equilibrium (KMS) state, is the Fock space 


(2.5) Hn = H^ = ^ 4'ym((R X F2)^ {du X d^r), 

n>0 


taken over the single-particle space L^(r x S‘^,du x dS), where dS is the uniform measure on 
pEU. carries a representation of the GGR algebra in which the Weyl operators are 
given by IF(E) = y^here 0(4) = ^(a*(4) -F a(4)). Here, a*(4) and a(4) denote 

creation and annihilation operators on J/g, smoothed out with the function 


( 2 . 6 ) 


Mn,S) = 


u 


\u 


1/2 


/(m, S), m > 0 

-7(-w,s), M<o 


1 - 

belonging to L^(r x du x dS). It is easy to see that the GGR are satished, namely, 
(2.7) Wif0)W{g0) = + g^). 


The vacuum vector G G J/g represents the inhnite-volume equilibrium state of the free Bose 
held, determined by the formula 

(2.8) {fi.ir(*)!J) = exp {-!(/, coth(/?|fc|/2)/)} . 

The Weyl algebra is represented on hF(/) t W{fp), for functions / G L^(R^) such 

that (/, coth(0|/c|/2)/) < cx). We denote the von Neumann algebra of the represented Weyl 
operators by W/g. 

The combined system-reservoir Hilbert space is then H, fl2.ll) . and the von Neumann algebra 
of observables is 


( 2 . 9 ) 


'm = B(Sjs)»ls®H>j CB(W). 
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The coupled dynamics is given by 

(2.10) a\A) = A eM. 

It is generated by the self-adjoint Liouville operator acting on Ti, 


(2.11) 

L 

= 

Lo +AJ 

(2.12) 

Lq 

= 

Ls + 

(2.13) 

I 

= 

V - JVJ. 

Here, Ls 

= Ps <H) Is — Is ® 77s and Hs is 

the 

system Hamiltonian acting on i^s- L^ = dr{u) is 


the second quantization of multiplication by the radial variable u. The interaction I in fl2.13l) 
is “in standard form”, involving a self-adjoint interaction operator V acting on "H and the 
modular conjugation J, which acts as 


(2.14) J{A 0 Is ® WiMu, S))) J = Is 0 A 0 iy(/^(-u, S)), 

where A is the matrix obtained from A by taking entry wise complex conjugation (matrices are 
represented in the eigenbasis of Hs). Note that by (12.61) . we have /^(—u, S) = fp{u, E). 

By the Tomita-Takesaki theorem [8], conjugation by J maps the von Neumann algebra of 
observables fl2.9p into its commutant. In particular, V and JVJ commute (strongly on a 
suitable domain). For more detail about this well-known setup we refer to [211 El El] and 
references therein. We have in mind two commonly used forms for V, 

(2.15) Vi = G 0 Is 0 0(h/3)h.c. or = G 0 Is 0 hF(h^)-h h.c., 

for some matrix G on i^s and where is a (represented) form factor, obtained from an 
h G L^(r^) by fl2.6l) . The interaction Vi is standard. V 2 comes about when considering the 
spin-boson system at arbitrary coupling strength [2l|, see Section [T2l 

The vector representing the uncoupled (a*,/3)-KMS state (A = 0) on is 

(2.16) flo.KMS = ^S,/3 ® 
where fls,/3 is given in fl2.4l) . 

In this setup of open systems, one can derive Condition (A2) from a global limiting absorption 
principle as follows. 

Theorem 2.1. Let Pr = ls0 |f2)(12|. Suppose that there is a dense set V G LL with Ran/Pj^ C 
V, s.t. for all ^V, 

(2.17) sup ^ G(0,i/;), /c = 0, ...,3. 

2GC- 

Then Condition (A2) holds with a = g/2 = mine^e'{|e — e'|}/2. 

We prove Theorem 12. II in Appendix lAl (The constant C{(j),'if) in fl2.17p may differ from that 
in (A2).) 
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2.2. The spin-boson model for arbitrary coupling strength. The spin-boson Hamilton¬ 
ian is [SB] 


(2.18) H = -|Acr^ -F \ea^ + H^ + \qocr^ ® 0(h), 

where ax and a^ are the Pauli matrices 






0 

-1 


and A, £ G R are the ‘tunneling matrix element’ and the ‘detuning parameter’, respectively. (We 
use units so that h takes the value one.) The reservoir Hamiltonian is i^R = \k\a*{k)a{k)d^k. 

The coupling constant is go ^ R- Associated to the Hamiltonian H is the Liouvillean 


(2.19) L — L^ + L-^ + qoI 
with 

(2.20) Ls = 1q2 — 1^2 (g) ihg, 

(2.21) Lr = dT{u) 

where Hg = —lAcr,j, -|- ^ea^, and / given in fl2.13p with V = ^a^ ® <p{h). The total Hilbert 
space is given by (12.ip . (12.3p with and (12.5p . 

In order to be able to analyze the spectrum of L for arbitrarily large couplings go G R, one 
applies the unitary (‘polaron’-) transformation U (see [21|) 

(2.22) U = uJuJ, where „ = 
resulting in a new Liouvillean [26] |24| 

(2.23) C = ULU* = £o + AJ, 
where 

(2.24) £o = = I (o'z (2) Is — Is + Tr 

and 


(2.25) / = -l(V-JVJ), V = a+®ls®W{2f0) + a-®ls®W{-2fp). 
Here, <7+ and cr_ are the raising and lowering operators and 

(2.26) f0 = {-^qoh/u)p. 


The non-interacting KMS state associated to Cq is 


(2.27) 


^0,KMS = <8 H, 


^S,/3 = 


Ve-/3A2 + Qf^e/2 

and the interacting KMS state associated to C is (Araki’s perturbation theory of KMS states) 

e-^(^“+^^)/2To,KMS 


(2.28) 


^KMS — 


|e-/^(^o+‘^^)/2To,KMs| 
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Note that the spectrum of Cq consists of a purely absolutely continuous part covering all of 
R, in which are embedded the eigenvalues e = ±£ (each simple) and the doubly degenerate 
eigenvalue e = 0. The following is the main result of [23]: 

Theorem 2.2 f[23]b Recall that (•)^ is defined in fl2.6p . Assume that (1 + \iduY'){\h/u )0 G 
L^(r X du X dS) for some rj > 2. Then, for any qo G R, go 7^ 0, there is a constant Aq > 0 
s.t. if0< |A| ^ Aq, then C has purely absolutely continuous spectrum covering R and a single, 
simple eigenvalue at zero. The associated eigenvector is Tkms; 112.28p . 


Remarks 1. The eigenvalues e = of Co are unstable, while e = 0 is partially stable. The 
associated (simple) eigenvalue of A is i? = 0. 

2. The spectral properties of L and C are the same, as the operators are unitarily equivalent 
to each other. 

3. The KMS state associated to L is given by 


(2.29) 


f^KMS — — 


^-l3(Lo+qo(T„mc2'S><t>(hp))/2Q^^^^^ 


where Lq = Ts(A) + Lr, see fl2.20p . fl2.2ip . One shows that Oo,kms is in the domain of 
g-/3(Lo+<?oo-zig)ic2i8)0(/i/3))/2 fQj, go, A G R (see e.g. [T^ 171 IS]!. 


We now verify assumptions (A1)-(A4) for the spin-boson system, i.e., for the operator C, 
{\2.2‘i\ . The eigenprojections of £o are given, for e G spec(£5),by 


Fe = l[As = e]®|0)(0|. 

Since l[£s = e](n‘± <8) ls)l[As = e] = 0 = l[£s = e](ls <8) cr±)l[£s = e], condition (Al) holds. 

To verify the limiting absorption principle (A2), let N = dr(l) be the number operator on 
Fock space (12.5p and put N = PfiN. Let A = dr{idu) and put A = Fr A. For a, z/ ^ 0, dehne 
the norms 


(2.30) II?IU. = I|A'‘'Ai + T)“''7I|. 

We have the following regularity properties of the resolvent 


Theorem 2.3. Let /i ^ 1 and suppose that dlfg G F^(r x du x dC), for j = 0, ..., 2/i -|- 1. 
We have 

(2-31) sup \£^{(j),R^^fi)\ -< \\(^)\\|,,2,^M\,,,2|^ 

z£C- 

(2-32) sup R^^fi)\ -< ||(/)|| 3 ,i||V'|| 3 ,i- 

2;GC- 


We give a proof of Theorem 12.31 in Appendix!^ The bound fl2.3ip with p = 4 implies fl2.17p . 
with the dense set 


(2.33) V = {fi : ||V^||4,8 < oo}. 
To see that Ran/Fa C F, it suffices to check that 

(2.34) IIA^(l + A^)W{2ff})Q\\ < oo. 
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It is not hard to use the relation 

(2.35) Idr(D), W{f)] = W{fm,Df) + i{/, Df)) 

(where (j) is the held operator, see the proof of Lemma 3.1 in [2l] for technical details) for 
D = —idu and D = 1 to see that fl2.34p holds provided G L^(r x S"^), j = 0,...,4. 
Therefore, Theorem 12.31 combined with this last observation shows that the assumptions of 
Theorem 12.11 are satished. Thus, by the latter theorem, assumption (A2) holds. 

Next, assumption (A3) is shown to hold in Theorem lA.ll flA.2p . Namely, the regularity in A 
of -^PelR^’^IPe is derived from that of -^PelR^IPe-, given in Theorem 12.31 fl2.32p . 

The Fermi Golden Rule Assumption (A4) is verihed by examining the level shift operators 
Aq and A±e. Aq is two-dimensional, given by (see [2l], Proposition 3.5) 

(2.36) Ao = 

where is the complement of Ps ,0 = |r2s,/3)(f^s,;3| in Ranl[£s = 0] and where oc 
(g) -g (g) ip_ (see (12.41) and (12.241) ). Also, 


(2.37) 

with 


poo 

/ dt cos(et) cos 


Jo 

TT 


Qi(^) = 


r-i = 


dbj —^sm(ci;t) and ( 42 ^) = 


du 


'0 


-^Q2(i) 


J(a;)(l — cos(a;f)) 


coth(/9a;/2). 


Here, J{uj) is the spectral density of the reservoir, dehned by 

(2.38) J{uj) = ^uj'^ ca > 0, 

the integral being taken over the angular part in R^. The function h is the form factor in the 
interaction fl2.18p FI Relation fl2.36p gives 


(2.39) 


Ao,o — 0 


and 


Ao,i = ir ^ 


Hence assumption (A4)(2) holds. The resonances A±e,o are the eigenvalues of the one-dimensional 
level shift operators A±^, which are easily calculated to be 


A±£,o = ±x + Air 


-1 


(2.40) 

where ±a; is the real part. Assumption (A4)(l) thus holds. 

Set (p+- = 99+ (8) etc. and 

^0 = w{fy)jw{fy)n, x* = w{f0yjw{fyrn, 

x+ = w{fyrjw{fy)n, x_ = jx+ = w{f0)jw{f0rn. 

The dynamics of the spin-boson system at arbitrary coupling is then explicitly given as follows. 


^The spectral density is related to the Fourier transform of the reservoir correlation function C{t) = 
Wfl,/3(e**^®y(h)e“®‘^^y(/i)) by J(a;) = ^ 7r/2 tanh(dw/2)[C'(a;) -I- C'(— w)]. Of course, it is assumed here, as 
it is in |26j , that the integral in (I2.37P does not vanish, so that r < 00 is a finite relaxation time. 
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Corollary 2.4 (Dynamics of the spin-boson system at arbitrary conpling strength). Suppose 
that ufj 3 ,difi 3 G L^(r X 3“^) for j = 0,...,4. For any Qq & R there is a Aq > 0 s.t. if 
0 < |A| < Aq then the following holds. 

Denote by fin = |DKMs)(f^KMsl the projection onto the coupled KMS state, fl2.29l) . Let 0, D G 
dom(LR) n dom(iVi^/2 ^ all t > 0 

(2.41) = (0,no0)+e‘*^^““(0,n'0) 

+ei*("+A"ae)^0^ n;^0) + n'_,0) + R{t), 

where R(t) -< | and 

oo = ir“^ 0(A), a±e = dux + + 0(A), 


(2.42) 

(2.43) n; 

(2.44) n'_. 


(efc/2 + 0 

-k-)U++l 0 a;)(x„| + 0 A;> Ao'l) + 0(A), 

|v+_)(v+_|0|x+)(;!r+| + o(A), 

|tJ_+)(tJ_+|®|Jf_)(X_|+0(A). 


Note that the dynamics in Corollary 12.41 is expressed with respect to the original Lionville 
operator L (not the unitarily equivalent £). 

Proof of Corollary \2.4\ We have (0, e'*'^0) = (f/0, Theorem 1 1.1 1 gives the resonance 

expansion for (f/0, e^*'^O0) provided f/0, f/0 G V s.t. CU(j),CU'ilj G V, where V is the set of 
vectors with finite || • ||4^8 norm (see after Theorem 12.31) .One can easily show the bound 

WUlfha^u -< ||0||2a,J/- 

We conclude from fl2.45p that if 0 G dom(A'^+^0 _|_ i)(^ 2 q: _|_ some z/ ^ 1/2 and a 0 0, 
and if {iduYfg ^ T^(r x S'^) for j = 0,..., 2q!, then ||f/0||2a,j/ < oo- 

We infer that if 0 G dom(LR) fl dom(A^^O _|_ 1)[A^ + 1), then ^ V and Clhf G V. (To 
see the latter inclusion, we use fl2.35p with D = u and proceed as in fl2.45p .i 
The operators IIq and are given by 

(2.45) n'o = U*Pi-^gPnU + 0(A), n0, = U*P±,U + 0(A), 

where P±s are the eigenprojections of Cq associated to the eigenvalues ±e, is defined after 
fl2.36p and Pr = Is <8 |D)(r2|. It is easy to calculate ® Ll = (p_|__ (8 W{fg)*JW{fg)JLl, 

which shows fl2.43p . Similarly one obtains (12.441) . For fl2.42p . we first calculate 

Ps,p = + e"^^^2)"^(e^^/V++)(‘/'++l “ \f++){F—\ “ \F—){f++\ + ■ 

Then we use (12.451) to arrive at (I2.42p . This completes the proof of Corollary 12.41 □ 


Being a KMS state, Dkms given in (12.291) is separating for 071, which means that OJI'Dkms is 
dense in R, where 071' is the commutant of 071, see [8]. Any (normal) state cu on 071 is given by 
a normalized vector 4/ G P via uj{A) = (T, AT). We introduce the dense set 

(2.46) Vo = dom(LR) n dom{N^'^/^ + 1)(A^ + 1). 
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The set of states cj arising from vectors in 

(2.47) {T G "H ; ||4/|| = 1, T = BUkms for some B G OJt', s.t. 5*4/ G 5o} 

is dense (in the norm of states on Tt). We call it the set of regular states, S^eg- We also 
introduce the regular observables, 

(2.48) 9ftreg = {a G 971 : A^kms ^ ^^o}- 
Let us denote the coupled equilibrium state by 

(2.49) a;KMs(^) = (f^KMS, ^f^KMs)- 

Corollary 2.5 (Return to equilibrium). For any Uq G S^^g, ^ £ ^Tlreg; t ^ 0, we have 

(2.50) |a;o(a*(^)) — <^kms(^)| ^ + 1/^ • 

The constant Ca,ljo depends on the initial state Uq and the observable A, but not on t,A. 

One readily verihes that all states of the form ojs <8 wr, where ois is arbitrary and oir is the 
reservoir equilibrium, belong to iSreg. Moreover, all observables on the system alone belong to 
OTlreg, since Orms e Vq. 


Proof of Corollary \2.A Let 47 = ROrmS) B G 971', be the vector representing uq. Since 
B commutes with a^{A) and Orms is in the kernel of L, we have q*'{A)'^ = Ba^{A)QKMS = 
Be^^^AQinMS- Thus, 

(2.51) a;o(a'(7l)) = (4/, e'*Me-'*^4/) = (5*4/, e'^MORMs)- 

Now we apply fl2.4ip and, using that flo = |Orms)(^kms|) obtain directly 02.501) . □ 


3. Proof of Theorem 11.11 
For z G C_, |Rez — e| ^ a we dehne the operator 
(3.1) R. = - PJR ^ pIPe . 

As, z ^ e and A —)■ 0, Re approaches the level shift operator Ag. More precisely, we have the 
following result, in which 5 is, recall, given by 01.181) . 

Lemma 3.1. There is a constant c such that 7/ |A|, |Rez — e| < cS, and z & C^, then 

1. All eigenvalues of A^ are distinct. Call them Oej = aej{z), j = 0,..., mg — 1. Each Oej 
satisfies \Xej — aej\ <5/2 for exactly one eigenvalue Xgj of Ae. 

2. The eigenvalues aej{z) of Az, and the associated Riesz projections Qj{z) are analytic in 
z G C_, iRe^; — e| < c5 and continuous as Imz —)■ 0_. They satisfy the bounds -^ae,j{z), 
■^Qj{z) -< 1 for fc = 0,..., 3, uniformly for |Rez — e| < c5 and Imz ^ 0. 

The simplicity of the spectrum implies the spectral representation 

me —1 

Az = ^ ^ aej{z)Qj{z). 
j=0 


(3.2) 
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Proof of Lemma [Ql When necessary, we display the A-dependence of by ^^(A). (Here, 
Ax for a: G R is understood as the limit of A^, as z ^ x, z G C_.) We have Ae(0) = and, by 
assumption (A3), 


(3.3) 


||Ae - A^(A)|| -< |A| + - e 


Assumption (A4) implies that 


(3.4) 


II (Ae 


C) 


- 1 | 


-< 


1 

dist(C, spec(Ae)) 


Using the standard Neumann series for resolvents, together with the estimates (13. 3 h and (I3.4p . 
yields 


(3.5) 


(A.(A)-C)-^ 


(Ae 


C)"' + o 


/ |A| + \z - e\ \ 

Vdist(C,spec(Ae))2y ’ 


provided |A|, 1^ — e| < cq dist(C,spec(Ae)), for some constant cq independent of A, 2;. Let Cj 
be the circle centered at Aej with radius 5/2 and dehne 

(3.6) 0,(A, z) = ^ / (T(A) - Q-^dC,. 

JCj 

Note that Qj(0,e) is the Riesz eigenprojection of Ae associated to the eigenvalue Aej-. Using 
fl3.5l) and fl3.6l) we obtain that \\Qj{A,z) — Qj{0,e)\\ < 1, provided |A|, \z — e\ < Cq5/2 and 
I A| + |z — e| < Ci5^/4, for some Ci independent of A, z. This proves point 1. of Lemma [3.II 
Next, 

(3.7) Q'=|.Q, = ^j((^(A)-0-‘Y(A(T(A-C)-'A 1 . 

Since af,j{z) = TtA^Qj^z) we get = Tr(A'^Qj(2;) + A^Q'jiz)) -< 1. The statements about 

the higher derivatives follow in the same manner. This shows point 2 and completes the proof 
of Lemma 13.11 □ 


3.1. Estimates for z in a vicinity of a partially stable eigenvalue e. We introduce the 
operators 


(3.8) Q, = -Q,(e)PJR^,. - Rt.oJPeQi(‘) + AR^„IP,Q,(e)PJR^-, 


-i0+ ’ 




where Qj{e) is determined by fl3.ll) and Relio+ limit of = (L — z) ^ as z approaches 

e through the lower half plane. Note that Qj{e) = Pgj + 0(|A|), see also fll.lUl) . 


Proposition 3.2. Let he the orthogonal projection associated to the eigenvalue Eg of L. 
We have, for 0, -0 G "P, 

me —1 

J^it) = e''®^(0,nE^0) + eihe+A2ae,,(e)) 

i=i 


(3.9) 


(*(*. {QAO + A^Qj)i)) + Reit) 
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where (recall that w > 0 is the arbitrary parameter in fll.l7p ) 

1 + e^yt 


\Re{t)\^C- 


t 


for a constant C independent of A,t,w, and where 

(3.10) aej{e) = Xe,j + 0(|A|), Q^e) = + 0(|A|). 

Proof of Proposition We apply the Feshbach map with projection Pg (having rank me), 

(3.11) ^(z) = ^(L - z;Pe) = e - z + A^A^, 

where Az is given in fl3.1l) . Due to Condition (A2), z Az is analytic for 2; G C_, \z — e\ < a, 
and its 2;-derivatives up to degree 3 stay bounded as Imz —)■ 0_. According to the decomposition 
fll.25p . we have 


(3.12) f e“‘{4>,iWi>)dz= j e>“ [(■(■, + + 

where 

(3.13) B{z) = -Adiz)-^PeIRy - AR^IPediz)-^ + A^RyiPe^{z)-^PeIRy. 

To examine 5(2;)“^, we use the spectral representation of the operator Az- 

3.1.1. The contribution to fl3.12p from ■ In this subsection, we will simply write 

Oj = Oej, 0 ^ J ^ TTle — 1, 

to ease the notation. Due to fl3.1ip and Lemma 13.11 


me —1 


(3.14) 


= E 


e — z 


+ A^a. 


-Qj- 


j=o - - • J 

We analyze the first term on the right side of fl3.12l) . using that 


(3.15) 


d 

dz 


me —1 

= E r, 

j=0 


where 


T- = 


1 - A^a' 


.z - e 


A^a, 


Qj + 




Q'r 


z + A^a,' ^ 


We examine the singularities of Tj in z. By the isospectrality of the Feshbach map, we know 
that e — Pe + ^^Ae^ has an eigenvalue zero (see also fll.SiB . Therefore, e — Pg + A^ao(Pe) = 0. 
Also, aj{z) = Aej + 0(1 A| + \z — e\) for j = 1,..., rUe — 1. Consider first Tq. We have 

2; - e - A'^ao{z) = z - e - A^ao(Pe) + A^(ao(Pe) - ao{z)) = z - Ee + A^(ao(Pe) - ao{z)) 


and so 
(3.16) 


1 - A^a'o _ _ 1 - A^a^z) _ 

{z-e- A2ao)2 {z - Pe)^ [1 - A2(ao(^) - ao{Ee))/{z - Pe)]^ 

1 _ 1 hjy 

(z-Pe)n-AX(^e) (Z-Ee)^' 
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where 

(3.17) 


h{z) 


_ 1 - _ 

[1 - A2(ao(^) - ao{Ee))/{z - Ee)Y 


1 - A^a(,(i^e) 

[l-AVo(Ee)]2 


0{^^\z-E,\^). 


To arrive at (I3.17p . we expand h{z) around a point zq G C_ which is very close to 


h{z) 


h{zo) + {z- zo)h'{zo) + 




ds'h''{s'). 


The integrals are over paths (straight lines) in the lower complex plane. Then, sending zq —)■ 
using that h{Ee) = h'{Ee) = 0 and controlling the double integral with the third derivative of 
h, we arrive at fl3.17p . In this argument, we assume the derivatives up to order three to have a 
continuous extension as Imz —?■ 0_. 

Thus 


(3.18) 


1 Qo{E,) 1 [Qo{z)-Qo{E,)]/{z-Ee) 

{z - E,y 1 - A2a'(Ee) z-E, 1 - A2a'(Ee) 

-i--+ o(A^). 

z- Eel- A2[ao(^) - ao{Ee)]/{z - Ee) 


An expansion of the sum of the second and third term on the right side of (I3.18p shows that 
this term is 0(1) uniformly in z G Oe, giving the bound 


(3.19) 
Therefore, 

(3.20) 



1 Qo{Ee) 

{z - EeY 1 - A2a' (Ee) 


e^*"(0,ToV^) 


(0,go(Ee)i/^) f _ 

1 - A^a'oiEe) Jg^ {z 


0(1 +A^). 

itz 

■3^<iz + 0(l). 


The remaining integral on the right side is estimated using the standard Cauchy formula from 
complex analysis. Namely, we complete Qe into a rectangular closed path, adding the vertical 
pieces C± = {eAa + iy : y G [—w, 72]} and the horizontal roof {x + iR : e — a ^ x ^ e + a}. 
Then we obtain from the Cauchy integral formula of basic complex analysis, upon taking 
7? —)■ oo, that 


(3.21) 


Atz 


dz = 27ri(e‘*^)'|2=£;^ + 0{e^^/t) = it27rie^*^'= + /t). 


Ig. (z - EeY 


The remainder term comes from the integrals along the two vertical pieces of the path, which 
are bounded above by J^e~'^*dy. Note that tc > 0 is arbitrary (see (11.161) 1 and we will take 
tc —?■ 0 which will make the remainder in fl3.2ip to be 0(l/t). Combining (13.201) and (13.211) 
yields 




(0,Qo(7^e)V^) 
1 - A^a'oiEe) 


+ 0(l/t + e“'7f2). 


(3.22) 


1 1 
it 27ri 
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Next, we analyze Tj in fl3.15l) . for j ^ 1. Recall that aj(e)|A=o = K,j are the eigenvalues 
with strictly positive imaginary part of the level shift operator Ag. The integrands behave in a 
different way now, since Imaj(e) > 0 (while before, ao(-E'e) = 0). The following bound is useful, 

\z — e — A‘^aj{z)\ = jz — e — — A^(aj(e) — Agj) — A^(aj(z) — aj(e))| 

^ jz — e — A^Aejl — ciA^(|A| + \z — e\) 

(3.23) ^ — e — A^Aejl, 


provided that |A|, \z — e\ ^ C2lmAej, where ci and C2 = l/(2ci) are independent of A and .2. 
To arrive at the last inequality, fl3.23p . we proceed as follows: the inequality is equivalent to 
2ciA^(|A| + 1^; — e|) ^ \z — e — A^Aejl- Now \z — e — A^Aejl ^ jlmz — A^ImAejl ^ A^ImAej, 


since Imz < 0 and ImAgj > 0. 


We have 
(3.24) 
where 


rj. ^ Qj(e) _ A^a'j(z)Qj(z) „ ^ 

' (z-e-A^aj(e))^ (z-e-A^aj(z))^ ^ 


z — e 


Qji^z') Qj(c) 


. _ _ 

^ z — e — A^aj{z) \_z — e — A‘^aj{z) z — e 


Q'iz) 


1 


z — e — A‘^aj(z) 
\z — el 


Qji^) - Qjie\ _ ^ A^aj(z) Qj(z) - Qj(e) 


z — e 


e — A^aj(z) z — e 


A2 


z — e — A^Aejl 


(3.25) = O ( 1+ 

and 

(3.26) 

^ -2A^(z - e)(aj(e) - aj(z)) + A^laj(e)^ - aj(z)^j^ ^ 

~ ^ A2^ . A2^ ^„M2 vRUj u 


/ A^jz — e 

P + A^l 


— e 

V \z — e 

-A^Aejl 

4 


[z — e — A'^aj{z)Y[z — e — A‘^aj{e)Y 

To arrive at the estimates fl3.25p and 03.261) we have used 03.23p . Similarly, the second term on 
the right side of 03.24p is 0{A'^/\z — e — A^AejP) and so we obtain 

Qj{e) 


(3.27) 
with 

(3.28) 


T, = 


{z — e — A‘^aj{e)Y 


+ Tj, 


mw -< 


z — e 


z — e 


A2A 


+ 


A2 


ejl 


Z — e 


A2A, 


+ 


A2| 


z — e 


TA^I 


z — e 


e,J I 


Z — e 


A2A 


ejl 


Note that, with Agj = ij + ir/j and w, r]j > 0, we have 


\z — e 


\z — e — A2Ae,j| 


^ 1 + 


A2|A, 


ejl 


Z — e — A2Aej| 


^ 1 + 


I A, 


ejl 




-< I5 


so 

(3.29) 


\z — e 


\z — e — A'^Aej 


-dx -< 1. 
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Also, 

(3.30) 


A2 


\z — e — A^A 


rdx < 2 


1 roo 

rdx ^ 2 


dy 


e,J\ 


x"^ + lS!^rd 


lo y + Vi 


-< 1 . 


Very similarly, one sees that 

(3.31) 

and it follows that 

(3.32) 


A2|^-e|2 +A^l^-el 


\z — e — A^A 


dx -< 1 


ejl 


dz^l. 


Next, using Cauchy’s integral formula as above ( (I3.19I) - (I3.21I) ). we obtain that 


Atz 


Qj{e) 


-dz = it27rie‘*("+^'“>("))g,(e) + 0{e^^/t). 


AT {z-e-A^a,{e)r 

Combining fl3.33p with fl3.27p and fl3.32p . we see that for all j > 1, 


(3.34) 


1 1 

it 27ri 


^^^A,TjAdz = + 0(l/t + e“'7t2). 


At this point, it is instructive to explain the coefficient (1 — A‘^aQ{Eg))~^ in front of the 
non-decaying term in fl3.22p . Let be the projection onto the embedded eigenvalue of L. 
We have (in the strong sense) LIe^ = lim£^o+(i£)(-^ — E^ + ie)“\ so by fll.251) . 


(3.35) 


Ue^ = hm (ie) {j?(L -E, + ie; P,)"' + B{E, - ie) + Pj.jJ . 

e^0+ ^ e ^ 


The Pe-block of the decomposition is, by fl3.14p . 
(3.36) 


me —1 


PellEePe = lim (ie)jJ(L - Pe + ie; Pe) ^ = hm ^ - 


le 


e — Pe + ie + A^aAPe — ie) 

j=0 •> 


Qj{Ee - ie) 


For j ^ 1, we have Imaj(Pe) > 0 and the corresponding term in the sum vanishes in the limit 
e —)■ 0. Hence 


(3.37) 


PeH^^Pe = lim (ie) 


go(Pe - ie) 


Qo{Pe 


0+^ ^ e — Pe-|-ie-|-A2ao(Pe — ie) 1 — A2ao(Pe) 


We have used the relation e —Pe +A^ao(Pe) = 0 (see after (I3.15p ). Therefore, the non-decaying, 
oscillating term on the right side in fl3.22p is 


(3.38) 


ZitEe 


(0,go(Pe)V’) = e‘‘®g0,PenE^PeA). 


1 - A2a'(Pe) 
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Finally, we combine fl3.22p . fl3.34p and fl3.38p to arrive at 

(3.39) +0(l/t + e"'7^^)- 

On the right side appears the Pg-block PJlE^Pe of the projection IIe^. The contribntions of 
the terms in (I3.12p with B{z) and add the remaining blocks to hnally give the full 

expression (0, 


3.1.2. The contribution to fl3.12p from From fl3.13l) and fl3.14p . we have 


(3.40) 

mg —1 

i{^,B{z)i,) = A ^ 


j=0 

where 


(3.41) 

q,{z) = [Q^PJR^p 


1 - A^a' 


■ (Ij + 


qp 


The summand in (I3.40p is of the same form as Tj in (I3.15p . with Qj replaced by qj. We may 
thus repeat the analysis leading to fl3.22p and (13.341) . giving 


1 1 


Ago(^e) 


me —1 


( 3 - 42 ) 77 ;^/ e“‘£{4',B{z)i,)dz = “TiPz- n + ^ 


it 27ri Jg^ ' 


1 - A2a'(Pe) 


i=i 


0{\A\/t+\Ap‘/r). 


Recalling fl3.35p . fl3.37p and using that R^l-ie stays bounded as e —)■ 0+, we get 


(0, = ^hm (ie)(0, {d{L - + ie; p.) ^ + B{E^ - ie)} 0) 


(3.43) 


= ( 0 , PellE.Pefj) + 


Ago(Pe) 

l-A2a'(Pe)' 


3.1.3. The contribution to fl3.12p from -^Rz‘'■ By assumption (A2), (11.Sp . we have (0, ^R^’^ip) -< 
1 and so we get 

(3.44) i2_ / 


it 27ri 


t 


Combining (I3.12p . (I3.39p . (I3.42p . (I3.43P and (I3.44p gives (13.9p . This concludes the proof of 
Proposition 13.21 □ 
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3.2. Estimates for z in a vicinity of an unstable eigenvalue e. The analysis is the 
same, actually somewhat easier, than the one presented in Section 13.1[ Indeed, for an unstable 
eigenvalue e, all the Xgj have strictly positive imaginary part (see Assumption (A4), (1)). 
Therefore, we can proceed as in Section [3Tl and in fl3.15p . all the terms Tj (even for j = 0) are 
now treated as above, after fl3.22l) . We immediately obtain the following result. 


Proposition 3.3. We have, for 0,-0 G 


me —1 


(3.45) 


j=0 

where (recall that w > 0 is the arbitrary parameter in (ll.lTh ) 


for a constant C independent of A,t,w and 
(3.46) 

Also Qj is defined in 


®ej(e) — Aej + 0(|A|), Qj{e) — PeJ + 0(|A|). 


3.3. Estimates for z away from the eigenvalues e. On the unbounded set Qoo we use the 
relation (ll.38p . This will help to ensure that the integrand is decaying sufhciently quickly at 
inhnity. Setting 0 = (L — i)0 and = (L + i)'0. 


(3.47) 

(3.48) 


Atz / 


{(j), j-R^Adz 



-2{cj>,R^A {<P,iRA) \ , 

{z + iy {z + iy J 



(z + l)^ {z + l)^j 


The integral fl3.48p is -< ||0|| \\L^p\\. The terms on the right side of fl3.47l) involving R^ and -^Rz 
are estimated using Assumption (A2)(2). Thus 


(3.49) 


1 1 
it 27ri 



j-^RzAdz -< 


t 


We have proven the following result. 


Proposition 3.4. We have 
(3.50) 


Joo{t) -< 


O((L-i)0, {L + i)fj) 
t 


3.4. Proof of Theorem II.11 We combine the estimates in Propositions 13.2113.31 and 13.41 to 
obtain the expansion fll.121) with the bound 


\R{t)\ -< 


1 + e^'Vt 

t 


for the remainder. Since ta > 0 is arbitrary, we have \R(t)\ -< 1/t. 


□ 
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Appendix A. 


A.l. From the global to the local limiting absorption principles. In this appendix, 
we derive the limiting absorption principle with projection P^, e G spec(Ls), from that with 
projection Pr given in Theorem 12.11 fl2.17p . Let 

= (PtLP,^ - rR^npi 

and recall the dehnition of the energy gap g, fll.lQp . 

Theorem A.l. Assume the conditions of Theorem \2.1[ Then, for 

(A.l) sup = C'(0,'0) < cx), /c = 0, ...,3. 

zSC_ : \z—e\^g/2 

Furthermore, fll.6p holds. Moreover, suppose that for G V, sup^g(;_ |^^((/), ^ 

(7(0, "0); k = 0,1. Then, for all e, 

(A.2) sup 1 ^( 0 , Pf ‘= 0)1 ^ ( 7 ( 0 , 0 ), fc = 0,l. 

2:GC_ 


(Note: the constants C{(j),'if) may differ from the one in (I2.17p .) 

Proof. We hrst show (lA.ip by expressing Pf® in terms of P^ and then using the bound 
fl2.17p . To do so, consider the operator 

(A.3) K = P^LPt + iP,, 

where i is the imaginary unit. We have {K — z)~^ = (i “ ^)~^p3 © Pf'" and therefore 
(A.4) Pf^ = {K - z)-^P,^ = P^K - z)-^P,^. 

Next, denote 


(A.5) dz = d{K-z;Pn). 

Then by (1B.2P (with (^ = Pr = Is © |f2R){fls|) 

(A.6) {K -z) ^ -ARP^P^IP^ 1 ) ( ^0 
where we have set 


1 -APJPj^pW 
0 1 


(A.7) Pe = l[As0e]PR = P,^PR. 

Due to (1A.7P and Theorem 12. ip V0 G P the maps ^ i-A PelP^Rf^fj and ^ i-A (0|Pf'^PR/Pe are 
three times differentiable with bounded derivatives for z G C_. 


Next, since Pe'^Pr = Pr and PrPe = 0, we have (P^PPj 




R 


I RanP, 


R 


-= Pf^ and so 


(A.8) 






AI - IR^^I) Pe. 


Since Pr/Pr -< 1 and Pb^IR^^IPr -< 1, we have 

(A.9) = (i - z)-^P, © (Ls - ^)-'Pe(l + 0(A)) 
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provided z E C- and [Rez — e| ^ ^g. The remainder term in flA.9p is uniform in these and 
thus, for A small enough, -< 1. We obtain -< 1 and, taking further 

z-derivatives, 

(A.IO) fc = 0,...,3. 

Combining (lA.4p . (1A.6p . (lA.lOp and the regularity (hrst three ^-derivatives bounded for z G C_) 
of the matrices to the left and right in flA.6D discussed above yields the result flA.lD . 

To prove fll.bp . the limiting absorption principle away from the eigenvalues of Lq, we apply 
the Feshbach map with projection Pr = Is 0 |r2R)(r2R|, 

(A.ll) = diL - z, Pr) = Pr(Ls + AI-Z- A^IR^^I)Pr. 

According to the decomposition fll.25p . 

(A.12) £R,i,) = {4,, id(z)-^4) + {4>, iB(z)i’) + (■#■. iR'Z-i’), 

where 

(A.13) B{z) = -A:S{z)-^PkIR'z^ - AR^IPn^S^z)-^ + A^R^^IPn^{z)-^P^IR^. 

For any 2; G Roo, we have ||(Ls — ^)~^|| ^ l/o; and therefore 

(A.14) = (Ls - [1 + APr(J - A/Pf-J)P r(Ls - P I- 

We then obtain at once from (lA.14p . (1A.13I1 and (12.171) that 

(A.15) |:5(z)-‘^l, {<l>,iB(z)i,)<i\^\C(<l>,i'), {4,iR^^)-<C{4,4)- 

The bounds flA.lSp together with flA.12p imply fll.bp . 

Finally, we prove flA.2p . We use again relations (IA.4P and flA.bp to express Pf® in terms of 
Pf®-. Taking the A-derivative in (IA.6P results in taking A-derivatives of Pf^ and The 
hrst one is controlled by the assumption in Theorem lA.ll the second one is controlled by flA.SD 
(as for the ^-derivatives above). □ 

A.2. Proof of Theorem 12.31 For ?7 > 0, we introduce the regularized Liouville operator (see 


(A. 16 ) 


C{r]) = Co + ir]N + AI{r]), with /(r;) = ( 27 r) / f(s)Tr,s{I)ds, 


and where Tt{X) = with A = dT{idu)- Here, / is a Schwartz function satisfying 

y(^)(0) = 1, fc = 0,1,... £(77) is a closed operator on dom(£o) Fl dom(A). 

The strategy of the proof is to derive estimates fl 2 . 3 ip . 02 . 321 ) for C replaced by the regularized 
£(77), namely, 

(A. 17 ) sup 

2:GC- 

(A. 18 ) sup 1 ^( 0 ,Pft^(77)1/^)I -< II^IIstIIV'IIs,!- 


2GC_ 


Here, -< means 01 .lip with a constant not depending on 77 > 0 . We have (</), R^^{ri)tjj) 
^(0, Pf'^V’) and ^(0, -E ^(0, Pf'^'0) as 77 —)■ 0+ (this follows from the fact that 


dzB- 
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Rziv) Rz{0) strongly as r; —)■ 0+, see [21], Lemma 4.3). Therefore, fl2.3ip . fl2.32p follow from 

dsiTi), dm. 

Throughout this proof, we will not indicate the dependence of C and I on rj and A. We will 
also simply write P instead of Pr. In particular, = {C{ri)—z)~^, where X = pixpi 

Let be an r^-dependent bounded operator in B{P'H) which is sufficiently regular in r; > 0. 
We dehne adyi(A^) = [A, and dX^ = — adA(A^). Note that OCq = 0 (recall that we 

write £o = R-oiji))- Since 5 is a derivation we have, according to Leibniz’ rule, 

(A.19) = -AR^{dI)R^, 

and for /i ^ 1, 

(A.20) Slflf]- = -A 

t = l 

Equations flA.lQp and flA.20p hold as equalities of bounded operators (note that dl is relatively 
A-bounded and rani?^ C dom(A), see [21], Lemma 4.3). It follows that 

(A.21) !;(<#', Iflfl'-iA) = {4>,!idAlR^]>‘ii) + {4>,a[R^ri,) 

(A.22) = {A4,, [R^ri,) - (|flf Ai,) + {4,, alBfl'-V). 

We are going to establish bounds on the right hand side. 

Proposition A. 2 ([21])- (1) Suppose the form factor f^ (see fl2.26p ) satisfies dffy G L^(r x 
S"^, du X dS), for j = 0, ...,£+ 1. Then, for i ^ 1, 

(A.23) \\N-R‘^{dI)N-R^\\ ^T]^. 

(2) VTe have 

(A.24) WnR^R^nR^W ^7]-\ 

(A.25) \\NR^R^ij\\ -< \\N-R\1 + A^YRfjW, 

(A.26) II^ + r^-^lliV(^-2)/V||, ^ ^ 2. 


In the arguments below in this proof, the biggest value of i in flA.23p we will use is £ = 2/i. 
Hence the regularity condition on fy in the Theorem 12.31 Combining flA.24p with flA.25p we 
get for j = 1, 2,... 


(A.27) 




Moreover, from (1A.26P we obtain for j,£ = 1,2,... 


P/2J 


(A.28) ||A^/2[^P]i^|| ^ ^-p/2j ||A^/2[^P]fc^|| ^ 


k=l 


k=l 


Since A ^ 1 and because of flA.27p we obtain for j, £ = 0,1,... 

(A.29) II A^/^[i?f]^'V'|| ^ r7-LA2J-i+i/2||jY-i/2('^ ^ AY^^i/jW + r/-LA2j || 
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From (1A.20I) and (1A.23I1 follows that 


(A^30) |(^,3|flf|'‘V)| A r,'5^ 

i=i 

Since (IA.27P holds with [i?f] replaced by [-Rf]*, we may apply flA.27p and flA.29p to flA.30p and 
obtain 

(A^31) \{^,d[R';]‘‘’P)\ ^ + 

J=1 

X ^„-[V2j-p+j-i/2|ljv-1/2(1 + A")‘''Ylt + ll/v'/Ylt 

By flA.22p we have 

(A.32) ^{4,, |flf l^iA) = (A4, - (|flf A4) + (4, d\R’’,Y4) 

and hence 

(A.33) \i^(4AR'!]'‘Y)\ < \\N-R^A4\\\\NRAR'AYn + \\NRAR'AY‘‘mN-R^m 


+ 11 *^ 111 , 2 ;i |I//'|Ii, 2 a(- 


By flA.27p we easily get 
(A.34) 


’1 

rj drj' 


- 11 + 1/2 




Since (0, [R^Yip) = (0, [R^ir] = 1)]^V') ~ [-Rf we conclude from flA.34p and 


from flA.27p for r; = 1 that 

(a.35) m[Km\ -< (i+r"+^/^)ii0iii,2.iiV'iii,2.. 

Let us now consider again flA.32p , but instead of using flA.27p for an upper bound apply flA.35p 
to (A0, and (0, [R^]^Aijj). In this way we get for flA.32p the upper bound 

(A.36) \i^{4, |iif|'‘UI ^ (1 + 'ri-'‘*^'AUh.2, m2.2,. 

Now we integrate (as before flA.35p i and obtain 

(A^37) imR/YY)] A (l+ir'‘+'''")l|.^ll2,2„ 11/^112,2„. 

The right side of (1A.37P has a power of rj reduced by one, as compared to the right side of 
flA.35p . but (j) and ijj contribute in a stronger norm. We continue this procedure to conclude 
the proof of flA.17p as well as the bound 


(A.38) 


sup 

zec- 


d dR- 


drj dz^ 


—{(j),R^^ili)\ -< ||0IU+i,2m ||'0IIm+i,2/2- 


We now prove the bound ( 1 A. 1811 . Note that 

(A.39) ^{4, R/4) = -{4, [R/lR/m- 
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We have (see flA.lQD ) 

(A.40) /i?f] = -A R^{dI)R^IR^ + R^{dI)R^ - A R^IR^{dI)R^. 

Next, we use flA.23p . (IA.24D and ||/|| ^ 1 to get 

(A.41) !(.#.,A .,'-'||JV'/2|flf]>||||flfV'll+r('l|iV''"[flfl>lll|iV‘'"flfV'll 
(A.42) 

By flA:^ we have 

(A.43) d[R^I A r^^-V^^^^^ll^lli.^llV'lliy- 

Recall that 


(A.44) = {A,/,, - ([flf/Bf]•,#., A^) + {,l,,dlR^IR^W. 

Due to flA.25p . the hrst and the second expression on the right side are bounded above by a 
constant times ?7“^||0||2,o||'0l|2,o- Thus we get from flA.43p . flA.44p 

(A.45) I A(0, [R^IR^mi A v-^<P\\2,i\m2,i- 

We integrate from 77 to 1 to obtain the estimate 

(A.46) A|ln(r7)|||0|h 


2,1 


2 , 1 - 


Next, we use flA.46p in flA.44p to get the better upper bound 


(A.47) 


I dr} 




3,1- 


Integration from 77 to 1 yields \{(j), [i?fIRf]!/))| -< ||0||3,i||'0||3,i, which, combined with flA.39p . 
implies flA.lSp . This completes the proof of Theorem 12.31 □ 


Appendix B. 

Let Q be an orthogonal projection and let Q-^ = 1 — Q. Let if be a densely dehned, closed 
operator satisfying RanQ C dom(iL). Let be the operator on RanQ-*-, given by Q^HQ^. 
For 2; ^ a{H^) we set R^ = — and we assume that ||i?^iLQ|| < cxo and HQiLRf || < cxo. 

Then we dehne 

(B.l) dz = d{H - z;Q) = Q{H - z- HQ^R^Q^H)Q. 

Theorem B.l (Feshbach-theorem, [5]). Let z ^ a{H^). 

(1) We have 

z e cr(iL) 0 e a(5'^) 
z e crp(iL) {0} 7^ ker5^ 

If z ^ (t{H) we have for Rz = (H — z)~^ that 

QRzQ = K\ Q^R,Q = -R^Q^HQ:Sf\ QRzQ^ =Q H R^,, 
Q^RzQ^ = R^+ R^Q^HQr,^QHQ^R^. 
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(2) If Q has finite rank, we have the additional characterization 

z G (t{H) ^ z E o'p{H) 4^ det(^5'^) = 0. 

(3) Let He the eigenspace of H for E G ap{H). The restricted projection Qlns ts a 
bijection from He to ker^^^;. Its inverse is 

0^0© -RfQ^HQfi. 


The previous theorem is proven in [S]. Point (1) of the theorem gives the block representation 


(B.2) 


Rz = 


-1 




0 


0 

r9 


-QHQ^Rf 

1 


1 0 

-R^Q^HQ 1 

Note that the three matrices on the last line are both invertible. Thus (recall fIB.ip ) 
1 QHQ^Rf 


H 


z = 


0 


:S{H-z-Q) 0 

0 Q^{H-z)Q^ 


1 0 

RfQ^HQ 1 


Remark. If Q' is another projection, satisfying Q'Q = Q' = QQ', then 
(B.3) d{d{H -z-Q)-Q') = d{H - z- Q'). 

Indeed, from point (1) in the theorem, we know that 

[dmn -z-Q)- g')]-' = - z- g)]-ig' = q'qr.qq' = Q'RzQ' = - z-, g')]"\ 

which implies fIB.Sp . 

Proposition B.2 (Weak Feshbach Theorem). Assume that H is self-adjoint, E G 'R\crp{H^), 
and that\mi^^Q^^{H — E + ie]Q) = — E;Q) exists (as a weak limit). If E is an eigenvalue 

of H with eigenvector fj, then 0 is an eigenvalue of '^{H — E;Q) with eigenvector Qfi. 

Proof. Applying the projections Q and Q-^ to the eigenvalue equation, we get for any e > 0 

(B.4) Q{H -E + ie)Qfi + QHQ^fi = ieQfi 

HQ'ijj + Q^{H — E -I ie)Q^'i/j = ie 


Applying the resolvent R-E-ie fo fhe second equality in flB.4D gives 
(B.5) Q-^fj = -RE-ie IIQ'i/j + ieR%_-^Q^fi. 

Using fIB.SD in the hrst equation of flB.4D . then letting e ^ 0_|_ and taking into account that 
RR%_ifQ^'f’ = E]Q^fi) = 0, gives — E;Q)fi = 0. (Note that Qfi 7 ^ 0 for 

otherwise, one easily obtains from the second equality in flB.4D that = 0 as well.) □ 
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